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Abstract 

Starting from the kinetic approach for a mixture of reacting gases whose particles 
interact through elastic scattering and a bimolecular reversible chemical reaction, 
the equations that govern the dynamics of the system are obtained by means of the 
relevant Boltzmann-like equation. Conservation laws are considered. Fluid dynamic 
approximations are used at the Euler level to obtain a close set of PDEs for six un- 
known macroscopic fields. The dispersion relation of the mixture of reacting gases is 
explicitly derived in the homogeneous equilibrium state. A set of ODE that governs 
the propagation of a plane travelling wave is obtained using the Galilei invariance. 
After numerical integration some solutions, including the well-known Maxwellian 
and the hard spheres cases, are found for various meaningful interaction laws. The 
main macroscopic observables for the gas mixture such as the drift velocity, tem- 
perature, total density, pressure and its chemical composition are shown. 
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1 Introduction 



The study of chemically reacting gases has been dealt with using the kinetic 
theory since the 1940s [1], and a large amount of literature has been developed 
about chemically reacting rarefied flows. This research field led to applications 
in the study of combustion [2], detonation (see [3] and references therein) and 
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atmospheric chemistry [4], fields in which a precise knowledge of the dynamics 
of the reacting gases mixture is required. 

Due to the increasing complexity that arises from a detailed microscopic de- 
scription, there are not many results in literature and they mainly concern 
numerical simulations and applications. 

Starting from the kinetic approach, fluid dynamic approximations have been 
constructed as simplified tools for a description, at a macroscopic level, of the 
mixture we are dealing with [5,6]. In order to recall some of the main con- 
tributions in this field, we can quote Ref.s [7,8] where the Chapman- Enskong 
approximation method was used in the homogeneous space case for an ac- 
curate evaluation of the chemical reaction rate and, more recently, the Grad 
method which was applied in both the homogeneous and non-homogeneous 
space cases, in order to calculate the reaction rate and transport coefficients 
[5,9]. 

In a recent paper [10], Rossani and Spiga obtained a more explicit derivation 
of kinetic equations with a chemical reaction, and succeeded in completely 
characterizing the equilibrium distributions and in proving an H-theorem for 
them. A generalization to many level molecules was performed in [11]. 

In the present work, the authors deal with a gas that undergoes to the following 
reversible bimolecular reaction 

A + B^C + D. (1.1) 

It should be mentioned that the kinetic theory governing reaction (1.1) has 
been the subject of investigation in recent years. For instance, in [12] the 
kinetic model of a diatomic gas with both dissociation and recombination 
reactions through a transition state was proposed, whereas in [13] the ex- 
tended Boltzmann equation governing reaction (1.1) was investigated in the 
Lorentzian gas limit. 

In the first part of this paper, starting from the nonlinear system of kinetic 
equations for a mixture of particles A, B, C, and D, we derive the balance 
equations for the different chemical species as well as the conservation equa- 
tions for the macroscopic quantities such as mass, momentum and total energy. 
By adopting the Euler approximation, in the case of a homogeneous equilib- 
rium state, we analytically derive the dispersion relation for the mixture of the 
reacting gases, where we found a relationship like k 2 P(iu) + uo 2 Q{iuo) = 0, 
where P and Q are polynomials of the first order in % u, with real coefficients. 
In the second part, we deal with the study of the propagation of a plane trav- 
elling wave. It can be remarked that, although many results can be found in 
literature about this important fluid dynamic problem for a dissociating gas 
(see [14,15] and references therein), the explicit treatment of a travelling wave 
with bimolecular reversible reactions is very recent [3]. 

In the dynamics of non-reacting gases, viscosity and heat conduction are re- 
sponsible for the structure of a travelling wave [15]. However, here we make the 
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assumption that relaxation due to elastic scattering is much quicker than the 
one due to chemical interactions. In this way we can neglect viscosity and heat 
conduction, so that only the chemical reaction is responsible for the structure 
of the travelling wave. Under the previous hypothesis on relaxation times, we 
construct a close set of nonlinear PDEs that govern the propagation of a plane 
steady travelling wave through the reacting gas. 

By taking advantage of the Galilei invariance of the Euler equations, we can 
adopt the reference frame in which the travelling wave is at rest. The resulting 
ODE system is reduced to a single nonlinear ODE which is solved numerically 
for several meaningful interaction laws, including the case of Maxwell and hard 
sphere molecules. Results on the whole gas (drift velocity, temperature, total 
density and pressure) as well as its chemical composition (molar fractions of 
the four components) are plotted. The effects of the different choices for the 
interaction laws on the structure of the travelling wave are discussed. 



2 Kinetic equations 



Let us consider a mixture of four gases A, B, C and D, which can interact, 
apart from all the elastic collisions, according to the reversible bimolecular 
reaction given in (1.1). In the sequel, we denote the particles A, B, C and 
D with the index % — 1, 2, 3, 4. The distribution functions fi(x, v, t) for 
particles i, endowed with masses m. and internal energy E i (available in the 
case of reactive interactions), obey the following system of nonlinear integro- 
differential Boltzmann equations (see [10] for details): 



^ + «.v/ t = m + QAf] , (2.i) 

where the vector / = / 2 , / 3 , f 4 ). In Eq. (2.1), Q z are the well known 
elastic collision integrals [15] while J i are the chemical collision terms given 
by: 



JAf](v) = J J 



"12 ' 34 
/^34 V 12 



dwdQ! , 

(2.2) 



where V 12 — \v — w\, V. u = \v' — w'\, ( — Q ■ Q', with f2 = (v — w)/V 12 , 
Q> = {v> - w')/V 34 and 



J. 



WW = // 



A*12 Ks4 



dwdfl' , 
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(2.3) 



where, now, V 34 — \v — w\, V 12 = \v' — w'\, ft — (v — w)/V 34 and Q' = 
(«' — w')/V 12 . The explicit expressions of the post-collision velocities v' and 
w' , (which can be found in Ref. [10]), are omitted here for the sake of brevity. 
The integrals J 2 [f] and J 4 [/] are obtained from J 1 [f] and J 3 [/] by introducing 
the following index exchange: 1 2, 3 ^ 4. In Eq.s (2.2) and (2.3), / = 
lll(V 34 , () and I** = I^(V 12 , () are the differential cross-sections of the for- 
ward and backward reaction, respectively, that are related by means of the fol- 
lowing microreversibility relationship: Vj (V 12 , () = n 3 2 4 Vj lH(V 3i , C), 
where /x r = m t m./M are the reduced masses, and M = m 1 + m 2 = m 3 + m 4 
(mass conservation) is the total mass of the reactants/products. The chemical 
collision terms satisfy the following important properties [10]: / J L [f]dv = 
J J 2 [/] dv = - J J 3 [f] dv = - J J 4 [/] dv = G34 - G12 , where 



Gu = fff V 34 ll 2 4 f 3 (v)f 4 (v')dvdv>dn> , (2.4) 
G 12 = HI V 12 f 12 /,(«)/,(«') dvdv'dQ' . (2.5) 

In the following we derive macroscopic balance laws starting from the kinetic 
equation (2.1). 

By integrating Eq. (2.1) over dv and using the relation / Q.[f]dv = [15], 
we obtain the balance equations for particles i: 



-Qf + ^-{n t u i ) = S i , (2.6) 

where the density n v the drift velocity u i of the species i, and the source 
term S t due to the chemical interactions, are given by: n i — J fX v ) dv, u. = 
(l/n.)fvf.(v)dv and S. = \ (G 34 - G 12 ) with \ = \ for i = 1, 2 and 
A. = -l fori = 3, 4. 

By multiplying Eq. (2.1) by m., m. v and Ei + m i v 2 /2, by summing over i 
and by integrating over dv, we obtain the conservation equations for the mass, 
momentum and total energy: 



§f + V.(„u) = 0, 

d 

— {pu) + V-{pu®u + lP) = 0, 

at 

de 

TTT + V 



dt 



{eM + lP)-u + q th + q h 



0, 



(2.7) 
(2.8) 
(2.9) 



where 1 is the unit tensor. Equations (2.7), (2.8), and (2.9) were derived by 
making use of the property [10] Yh=i I V , i ( t; )(^[/] + QAf]) dv = which is 
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valid for ipi = rrii, — m i v (three components), and ip i = E i +m t v 2 /2 (mass, 
momentum, and total energy conservation). In Eq.s (2.7)-(2.9) p = J2i m i n i is 
the total density and u = (1/p) ^Z i m.n i u. is the drift velocity of the mixture. 
We also posed IP = J2i m i I( v — u ) ® ( v — u ) fX v ) dv, q = \ J2i m l I( v ~ 
u) 2 (v - u) f.(y) dv and q. nt = Ei E. f(v - u) f.(v) dv = Ei E. n t («. - u) 
which represent the stress tensor, the thermal energy flux and the internal 
energy flux, respectively. Finally, e = \ p u 2 + e th + e int is the total energy den- 
sity where e th = | n x T = \ tr IP and e. nt = J2i E i n i are the thermal energy 
density (n = J2i n i is the total density, x is the Boltzmann constant and T is 
the absolute temperature) and the internal energy density, respectively. 
Equations (2.6), (2.8) and (2.9) are exact but do not constitute a closed sys- 
tem. After suitable assumptions, this system can be made close in order to 
obtain a simplified tool for the macroscopic description of the mixture, in 
terms of n t , u and T. 



3 Euler approximation 



If the characteristic time of relaxation due to elastic interactions is much 
shorter than the one due to the reactions, we can set Q i [f] = 0, whose solutions 
are the Maxwell distribution functions [15]: 



fAv) =n. 



m. 



3/2 



exp 



m. 



2*T 



(v - uf 



(3.1) 



At this point, u x , IP, q h and q. nt , as well as the source terms S v can be 
calculated as follows: 



where p = n x T is the pressure and 



(3.2) 



S. = A. 



z/ 34 (T) n 3 n A - u 12 (T) n, n 2 



where the effective collision frequencies are given by 



(3.3) 



X 



/ y /m 1 m 2 
J J Jv i2 f 4 2 (V 12 ,() exp 



m 1 v 2 + m 2 w 2 ^ 



2*T 



dvdwdQ,' , (3.4) 
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r i2 /T , ( m.,v 2 + m.w 2S 



V 3 J 34 (V 34 , C) exp ^ 3 2xT 4 I cfodtodO' , (3.5) 



and satisfy the following relationship [16]: 



^= K exp(^), (3.6) 
^(T) \XTJ 

(the Arrhenius law), where TZ = (/ii 2 //i34) 3/2 and AE = E 3 + E A - (E x + £ 2 ) 
(without any loss of generality, we can take AE > 0). At chemical equilibrium, 
G 12 = G 34 , that is, S i = 0. 

In this work we will consider the following class of microscopic collision fre- 
quencies: V 34 I 34 (V 34 , () = C 34 (()V.J where C 34 is a certain function of the 
cos(C) and 8 < 1 describes the interaction law. The cases 5 = and 6 = 1 
are well known as Maxwell and hard sphere models, respectively. After some 
calculations, from Eq. (3.5), we explicitly obtain 



^34 

= A 34 { X Tf\ (3.7) 

where Y is the gamma function [17], with 6 > —3. 
Due to the Arrhenius law, we soon obtain 



MT) = ^(xTpexp(-^) , (3.8) 
where A 34 is a constant. 



4 Dispersion relation 

Let us introduce the following notations to indicate a homogeneous equilibrium 
state n i = n*, T = T c with S i = = and we assume, without any loss of 
generality, u = 0. 

We consider the propagation of a harmonic plane wave: 

n i = n° + M i exp[? (x x ~~ u t)] , (4.1) 
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u = U exp[i (kx — cut)] , 

T = T C + T exp[i (k x - out)} , 



(4.2) 
(4.3) 



where u — u ■ i, where i is the unit vector along the wave propagation. 
By inserting Eq.s (4.1)-(4.3) into Eq.s (2.6), (2.8) and (2.9), after the lineariza- 
tion with respect to H i , U and T, we obtain: 



-i ou J\f x + % k n\ U = S\ , 
-iuu (A/; -N 2 )+iK (nl-nl) U = 0, 
-iuo (JVj +Af 3 )+iK (nl + nl) U = 0, 
-icu (jV; +Af 4 )+iK (nl+nl) U = 0, 

-ioup e U + ik X (n c T + T c =0 , 

+iK(^ X n e T e + ^E i n^U = 0, 



—iou 



where 



S* = u 34 (T c 



i/ (T c )< nj-i/,(r)n> 

34 V / 3 4 12 V / 1 



(4.4) 
(4.5) 
(4.6) 
(4.7) 

(4.8) 



(4.9) 



(4.10) 



Since at equilibrium u 34 (T°) n e 3 n e 4 = u 12 (T e ) n e 2 = R, from Eq. (3.6) it follows 
that *4(T C ) n e 3 n e 4 - v' 12 {T e ) n\n\ = —RAE/xT e2 . 

Non trivial solutions of the system of Eq.s (4.4)-(4.9) for J\f v U and T, are 
obtained if, and only if, the determinant of the coefficients vanishes. This 
condition gives the dispersion relationship, which links k and ou [16]: 



UU , 2 / OU 
T + k —%0UT\ = 



k 



0, 



(4.11) 



where the coefficients a c and at are the equilibrium flow and frozen flow speed 
of sound, respectively, while r is a characteristic time of the reaction. Explic- 
itly, we have 



5 n e T e 



% = 77 X 



Y 3n c VX^V 



(4.12) 
(4.13) 



7 



i = « 

T 



1 2 / AE \ 2 



(4.14) 



It is easy to verify that a,{ > a e , according to non-equilibrium thermodynamics 
[18]. We can observe that the dispersion relation can be written as k 2 P{i u>) + 
a; 2 Q{iuj) = where P and Q are first order polynomials of iuj. 



5 Travelling wave solutions 



When the macroscopic fields n v u i and T only depend on one spatial coordi- 
nate x and on the time t the macroscopic balance laws in the Euler approxi- 
mation become: 



at ox 

-^- t (pu) + ^-(pu 2 + n9) = 0, (5.2) 

~/t + dx' [ie + nd)u] = > (5 ' 3) 
where 6 = xT . 

In the following we consider the travelling wave solutions of the system de- 
scribed by Eq.s (5.1)-(5.3), that is, we look for solutions that only depend on 
the independent variable £ = x — vt where, due to the Galilei invariance of 
the system, we can set the velocity v of the travelling wave equal to zero. Eq.s 
(5.1)-(5.3) become: 



^ K «)=S i , (5.4) 
^(pu* + n6) = 0, (5.5) 
±[( e + n6)u} = , (5.6) 

which is a system of six nonlinear ODEs that have to be solved with the 
appropriate initial conditions. In the following we use a plus/minus superscript 
to indicate the asymptotic values at £ = ±oo. We anticipate that the soliton 
solution will connect two asymptotic equilibrium states. This means 



Sf = . (5.7) 
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We fix the values of the density n. and the drift velocity u as initial data, 
at £ — > — oo. Then, from Eq. (5.7), using Eq.s (3.3) and (3.6), it follows that 



9 = AE < log 



, n n 

n- 1 



™3 U 4 



(5.8) 



Taking into account the definition of S v Eq.s (5.4) can be rewritten as follows: 



d 
d 



[n 1 u) = S 1 , 



3 = 2, 3, 4 . 



(5.9) 
(5.10) 



Eq.s (5.5), (5.6) and (5.10) are of immediate integration and give the constants 
of motion of the system. Eq. (5.9) is the only dynamical equation and allows 
us to obtain the structure of the wave after integration. 
From Eq.s (5.10) we obtain: 



[n 1 (0-\n j (0]u(0 = (n;-\ j nj) 



u 



(5.11) 



and taking into account the definitions of the quantities n, p, e int , given in 
section 2, we obtain the following relations: 



n(0 = 

p(0 = 



«(0 ' 

<p 



«(0 ' 

e int (0 = ^y-A J En 1 (0, 

where u = n~ u~, ip = p~ u~ and e = e~ u~ + AE u~ n~ . 
By using Eq. (5.5) we obtain 



(5.12) 
(5.13) 
(5.14) 



*(0 = — [t-¥>«(0], 



(5.15) 



with t = p~ (m~) 2 + n _ which gives the temperature ^(£) as a function of 
velocity «(£)• 

If we now introduce Eq. (5.15) into Eq. (5.6), and use Eq.s (5.12)-(5.14) we 
obtain, after integration, 



a. 



«(0 



i = l, 



(5.16) 
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(it should be observed that rii is singular for u — 0), with a t = X t [(e— e)/AE — 
\n- - \ n-)u-\, (3 = -2ip/AE, 7 = 5r/2A£ where e = [p-(u-) 2 /2 + 
5n~ 8~/2 + e~J u~ . Eq. (5.16) gives the concentrations n^O as functions of 
only the velocity u(£). 

Finally, by using Eq. (5.16) in Eq. (5.9) follows the nonlinear ODE: 



du 



R[u\ 



with 



(5.17) 



R[ u ] = o f lN . (5.18) 

Eq. (5.17), after numerical integration, gives us u(£) and, if Eq.s (5.15)-(5.16) 
are used, we can find and (£). 

Some conditions are now discussed that have to be satisfied to obtain physi- 
cally meaningful results. It is well known that solitons of Eq. (5.17) reach finite 
asymptotic values which are two consecutive zeros of R[u] [19]. This implies, 
as we anticipated, S = 0, so that the solitary wave connects two asymptotic 
states where the system is at chemical equilibrium. By choosing the initial 
data n~ > and u~, we obtain 9 from Eq. (5.8), which must be positive. 
This implies the following constraint on n~ (the first condition): 



n n , 
-^^>TZ. 5.19 

n; n- 

Subsequently, we solve the transcendental equation S^u] = 0. Its roots are 
the asymptotic values ■u ± of u(£). Regular solutions of Eq. (5.17) require that 
there are no singularities of R\u] in the interval \u . , u m ], where u . = 

o L J L min* max J* mm 

mm(u~, u + ) and w max = max(«",« + ). Singularities in R[u] arise both for u = 
and 2 (3 u + 7 = 0. The condition u 7^ implies u~ u + > (the second 
condition), while 2/3« + 7 7^ for m G [ M min ) M nai ] is the third condition 
that has to be satisfied. Finally, by using Eq.s (5.15) and (5.16), from u + we 
obtain the asymptotic values n+ and 9 + , which must be positive (the fourth 
condition). When all these conditions are satisfied, by choosing an arbitrary 
initial value u G [u min , w max ], we can finally forward and backward integrate 
Eq. (5.17) in order to obtain the structure of the travelling wave which connects 
the two chemical equilibrium states at ±00. 

In the following we present and discuss some numerical solutions which connect 
the same asymptotic equilibrium, for various choices of 5. The behaviour of 
the gas is depicted in figure 1 as a whole, by plotting: the drift velocity u, 
temperature 9, total density n and pressure p. It should be observed that, as 



10 



-1,5 0,0 1,5 3,0 -1,5 0,0 1,5 3,0 




Fig. 1. Plot, in arbitrary units, of the drift velocity u, temperature 6, total density n 
and pressure p, for different values of the parameter 5: 5 = 1 (dash-dot line), 6 = 
(dash-dot-dot line), 5 = — 1 (dot line), 5 = —2 (dash line) and 5^—3 (solid line). 

6 < 1 in our example, for lower S is the higher are u 12 and z/ 34 , and the steeper 
the curves. On the other hand, in the case 6 > 1, we would get smoother and 
smoother curves for increasingly lower values of 5. 

The molar fractions x t are plotted in figure 2 for the four species, which are 
defined by x i = njn. This quantity has another interesting meaning. It is 
easy to realize that Xi is proportional to the current density j i = n i u. The 
slopes x[ = dXi/dt, are simply related by x[ = x' 2 = ~x' 3 = ~x' 4 , for each 
value of 5. 

Finally, it is possible to observe that, as the temperature decreases, z/ 12 de- 
creases with respects to u 34 . This means that the backward reaction prevails 
over the forward one, that is, the number of particles 1 and 2 increases while 
the number of particles 3 and 4 decreases, as can be seen in the plots. 



6 Conclusions 

We have studied a kinetic model for a mixture of gases with bimolecular re- 
versible reactions, in the extended kinetic theory. The kinetic equations have 
been obtained starting from the Boltzmann equation and taking into account 
both the elastic collision integrals and the chemical collision terms between 
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-1.5 0,0 1.5 3.0 -1,5 0,0 1,5 3,0 

% % 

Fig. 2. Plot, in arbitrary units, of the molar fractions % i5 for different values of the 
parameter 5: 5 = 1 (dash-dot line), 6 = (dash-dot-dot line), 5 = — 1 (dot line), 
5 = — 2 (dash line) and 5 — > — 3 (solid line). 

the particles. A closed set of equations is obtained in the Euler approxima- 
tion, under the assumption that the relaxation due to the elastic scattering is 
much quicker than the one due to the chemical interaction. Six partial differ- 
ential equations are derived which govern the evolution of the corresponding 
unknown macroscopic fields (the densities n t , momentum and total energy). 
Some numerical results have also been presented for the simple case of a plane 
steady travelling wave solution, for various meaningful interaction laws. Fi- 
nally, the shape of the soliton that represent the main macroscopic observables 
for the gas mixture (the drift velocity, temperature, total density, pressure and 
chemical composition) are discussed. 
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